Pinch Resonances in a Radio Frequency Driven SQUID Ring-Resonator System 
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In this paper we present experimental data on the frequency domain response of a SQUID ring 
(a Josephson weak link enclosed by a thick superconducting ring) coupled to a radio frequency (rf) 
tank circuit resonator. We show that with the ring weakly hysteretic the resonance lineshape of 
this coupled system can display opposed fold bifurcations that appear to touch (pinch off). We 
demonstrate that for appropriate circuit parameters these pinch off lineshapes exist as solutions of 
the non-linear equations of motion for the system. 

PACS numbers: 05.45.-a 47.20.Ky 85.25.Dq 



INTRODUCTION 



It is well known that the system comprising a SQUID 
ring (i.e. a single Josephson weak link enclosed by a 
thick superconducting ring) , inductively coupled to a res- 
onant circuit [typically a parallel LC radio frequency (rf) 
tank circuit], can display very interesting dynamical be- 
haviour P, 0, IE S Hi ■ This behaviour is due to the non- 
linear dependence of the ring screening supercurrent (/g) 
on the applied magnetic flux ($x) which originates ulti- 
mately from the cosine (Josephson) term in the SQUID 
ring Hamiltonian 0. In the past, limitations in the noise 
level, bandwidth and dynamic range of rf receivers used 
to investigate this coupled system have restricted the 
range of non-linear phenomena that could be observed, 
i.e. only relatively weak non-linearities in Is ($x) could 
be probed. However, recent improvements in rf receiver 
design have changed this situation quite radically. 

The regimes of behaviour of SQUID rings P| are 
usually characterized in terms of the parameter /? = 
(27rA/c) /<&o, where 1^ is the critical current of the weak 
link in the ring, A is the ring inductance and $0 {— h/2e) 
is the superconducting flux quantum. Thus, when /3 < 1 
(the inductive, dissipationless regime) Ig is always single 
valued in <I>x. By contrast, for (3 > 1 (the hysteretic, dis- 
sipative regime) this current is multi- valued in $x- Both 
these regimes can lead to very strong non-linearities in 
Is (^'x) lEQ' the inductive regime these will be 
strongest when (3 —>■ 1 from below, at which point Is 
becomes almost sawtooth in $x (modulo $0 ). For the 
hysteretic regime the existence of bifurcation points in 
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Is, at particular values of $x, can lead to even stronger 
non-linearities. In the absence of external noise these 
bifurcations, linking adjacent branches in Is ($x ), occur 
when $x — ±A/c. 

Although the behaviour of single weak link SQUID 
rings (as shorted turns) can be monitored at any non- 
zero frequency, most studies have been carried out at ra- 
dio frequency (rf, typically w 20MHz) for which a wide 
range of high performance electronics is available. We 
have adopted this frequency range in the work described 
here. In practice, signal to noise is improved (as in this 
work) by driving the SQUID ring through the intermedi- 
ary of a resonant circuit, almost invariably a parallel LC 
(tank) circuit with a high quality factor (Q). A schematic 
of the SQUID ring, inductively coupled to an rf tank cir- 
cuit (with parallel inductance and capacitance Lt and Ct, 
respectively), is shown in figure 1. In the most commonly 
used experimental scheme the dynamical behaviour of 
the ring is followed in the time domain through the so- 
termed SQUID magnetometer characteristics. In these, 
the rf voltage (Vout) developed across the tank circuit is 
plotted as a function of the rf current (/in) used to excite 
it, where /m is linearly amplitude modulated in time (fig- 
ure 1). This current generates an rf flux (peak amplitude 
ip) in the tank circuit inductor (coil) , a fraction of which 
{fi = M/Lt for a coil-ring mutual inductance of M) then 
couples to the SQUID ring. The supercurrent response 
in the ring then back-reacts on the tank circuit, and so 
on. Any additional static or quasi-static magnetic flux 
((f>xstat) that is required is usually supplied via a second 
coil coupled to the SQUID ring, again as depicted in fig- 
ure 1. For small /3- values (i.e. from < 1 to a few) the 
V^ut versus /in dynamics contain as their main element 
a set of voltage features repeated periodically in /in at 
intervals oc $o/A. The precise form of the SQUID mag- 
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FIG. 1: Schematic of the inductively coupled SQUID ring- 
tank circuit system. 



netometer characteristics depends both on the type of rf 
detection employed (e.g. phase sensitive or diode) and 
the value of $xstat- 

The non-linear screening current response can also be 
made manifest in the frequency domain. This can be ac- 
complished at a fixed level of (peak amplitude) rf flux 
jjnp at the SQUID ring, without the need for amplitude 
modulation. For example, in the inductive regime, with 
(3 very close to one from below, we have demonstrated 
that, at sufficiently large values of fold bifurcations 
develop in the resonance lineshape of a SQUID ring-tank 
circuit system. Indeed, using a bidirectional sweep (up 
and down in frequency) , and at certain values of pLip and 
^xstat, the system can even display opposed (hammer- 
head) fold bifurcations in this lineshape, one on each 
side of the resonance 0- Opposed bifurcations are not 
only most unusual but, to our knowledge, the possibility 
of their existence is scarcely mentioned in the literature 
on non-linear dynamics Q . In the case of the inductive 
SQUID ring, with (3 very close to one, we have found 
that these bifurcations are always well separated in fre- 
quency. However, in this work we will show that in the 
weakly hysteretic regime of SQUID behaviour (/3 w few) 
this separation can become extremely, perhaps vanish- 
ingly, small. Throughout this work we shall term this a 
pinch resonance. 

The one sided (above or below resonance) weak fold bi- 
furcation 9] is taken as the text book example of the on- 
set of non-linear behaviour in a strongly driven oscillator. 
The development of opposed fold bifurcations (on both 
sides of the resonance) in any non-linear system is an- 
other matter altogether. As far as we are aware, opposed 
fold bifurcations had not be seen experimentally prior to 
our report of hammerhead resonances in the frequency 
domain of coupled inductive SQUID ring-tank circuit sys- 
tems j7j] , and only once in the theoretical literature '^l . In 



itself, we consider that this constitutes a significant con- 
tribution to the body of knowledge concerning non-linear 
dynamical systems. That the cosine non-linearity of the 
SQUID ring can also generate, in part of its parameter 
space, something as exotic a a pinch resonance is simply 
remarkable, particularly since such resonances appear as 
solutions to the established equations of motion for the 
ring-tank system in the weakly hysteretic regime. What 
this result, and others, point to is that the non-linear dy- 
namics of SQUID rings coupled to other, linear, circuits 
is likely to prove very rich indeed [Tol |. the more so since 
the parameter space that has so far been explored is still 
very limited. As such, these coupled systems are likely 
to continue to be a focus of attention by the non-linear 
community, both in theory and experiment. 

With the SQUID ring treated, as in this paper, quasi- 
classically, there would appear to be many potential ap- 
plications of its intrinsic non-linear nature when cou- 
pled to external circuits, for example, in logic and 
memory (single ,11.] or multi-level [13 )i stochastic res- 
onance |l3j| and in safe communication systems based 
on chaotic transmission/reception However, in the 
last few years there has been a burgeoning interest in 
exploring the properties of SQUID rings in the quan- 
tum regime for use in quantum technologies, particu- 
larly for the pur pos es of quantum computing and quan- 
tum encryption |l5l IT^ ItR IT^ . This interest has been 
boosted by recent experimental work on quantum super- 
position of states in weak link circuits 0, |23, |^ 
and very strongly in the last year by reports in the lit- 
erature of experimental investigations of superposition 
states in SQUID rings 1^ ■ Quantum technologies 
will, of necessity, involve time dependent superpositions 
for their operation. Where quantum circuits such as 
SQUID rings are involved, this means the interaction 
with external electromagnetic (em) fields or oscillator cir- 
cuit modes. Correspondingly, the effect of applied em 
fields (classical or quantum) of high enough frequency 
and amplitude is to generate quantum transition regions 
where energy is exchanged between the ring energy lev- 
els. Given the cosine term in the SQUID ring Hamilto- 
nian _6J , these energy exchanges in general involve (non- 
perurbatively) coherent multiphoton absorption or emis- 
sion processes . These exchange regions extend 
over very small ranges in the static magnetic flux <&xstat 
applied to the ring (typically over 10"'^ to 10~''$o), com- 
parable to the width of the anticrossing regions reported 
in the recent exp eriments on superposition of states in 
SQUID rings 23, 24]. Such narrow flux widths imply the 
concomitant existence of extremely strong non-linearities 
in the screening supercurrents flowing in the ring. In gen- 
eral these current non-linearities generate non-linear dy- 
namical behaviour in external (classical) circuits used to 
probe the superposition state of the rin g. A s we have em- 
phasized in previous publications [23, IS US , following 
these non-linear dynamics is one way of extracting infor- 
mation concerning the quantum state of a SQUID ring, 
or any other quantum circuit based of weak link devices. 
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From our perspective, therefore, it is of clear importance 
to the future development of quantum technologies based 
on SQUID rings and related devices that these non-linear 
dynamics be studied in great detail. Given the results 
presented in this paper, the non-linear dynamics of cou- 
pled SQUID ring-probe circuit systems can display rich 
and yet unanticipated behaviour. With these rings oper- 
ating in the quantum regime this is even more likely to 
be the case. 



II. EXPERIMENTAL PINCH RESONANCES 

Although sometimes unappreciated, probing the non- 
linear dynamical behaviour of SQUID ring-tank circuit 
systems makes strong demands on the high frequency 
electronics used. In order to open up new non-linear 
phenomena in these systems we required rf electronics 
which combined ultra low noise performance, large dy- 
namic range and high slew rate capability. In practice 
we achieved this by using a liquid helium cooled (4.2 
Kelvin), GaAsFET-based, rf amplifier (noise tempera- 
ture ^ 10 Kelvin, gain close to 20dB) followed by a state 
of the art, low noise, room temperature rf receiver. With 
these electronics we could choose diode or phase sensitive 
detection of V^ut- Since phase sensitive detection gener- 
ally provides better signal to noise, this was always our 
chosen mode of operation for acquiring SQUID magne- 
tometer characteristics. 

In the work reported in this paper the SQUID rings 
used were of the Zimmerman two hole type [2I|. These 
were fabricated in niobium with a niobium weak link 
of the point contact kind. This weak link was formed 
by mechanical adjustment between a point and a post 
screw, the latter having a plane machined end face which 
had been preoxidized prior to its insertion in the SQUID 
block. The weak link was made in situ at 4.2 Kelvin 
via a room temperature adjustment mechanism. As we 
have shown in previous publications JJ^, the control we 
can exert using this mechanism is very good, sufficient 
in fact to make reproducible point contact SQUID rings 
with any desired (3 value from inductive to highly hys- 
teretic 

In figure 2 we show what we term the equal amplitude 
magnetometer characteristics for a weakly hysteretic (/? 
close to 2) niobium point contact SQUID ring coupled to 
an rf tank circuit. In these characteristics, taken at 4.2 
Kelvin with a static bias flux of n$o, the frequency of the 
drive current was set to the coupled SQUID block-tank 
circuit resonant frequency before a weak link contact was 
made (i.e. Ci;in/27r) and the voltage components {V^^^^ and 
V^^f. ) were maintained exactly 90° out of phase. How- 
ever, their relative phases with respect to Jin were ro- 
tated electronically until the average slopes of the two 
characteristics are essentially equal. This approach was 
adopted for convenience - it presented the periodic fea- 
tures in /in of the orthogonal phase components of I^ut 
at essentially the same amplitude. Here, the (measured) 
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FIG. 2; Experimental equal amplitude (phase sensitive de- 
tected and 90° apart in phase) rf SQUID magnetometer char- 
acteristics {Vout versus /in) of a weakly hysteretic (/3 = 2.26) 
niobium point contact SQUID ring at a static flux bias of 
$^,t^t = n^o- Here, A = 6 x 10'^° H, ujin /2n = 20.122M Hz, 
Lt = 6.3 X 10"*//, Q = 1058, K'^ = 4 x 10"^ and T = 4.2 
Kelvin. 



ring inductance A = 6 x 10 ^'^H, the measured ring- 
tank circuit coupling coefficient ^= = 0.004, 

Lt = 6.3 X 10"^//, Q = 1058 and the drive (excitation) 
frequency for /in, Win/27r = 20.122MHz. It is apparent 
from figure 2 that even in this weakly hysteretic regime 
the observed SQUID magnetometer characteristics are 
highly non-linear and contain a wealth of detail. 

Some of the time domain information of figure 2 can be 
recast in a quite dramatic form in the frequency domain. 
To extract this experimentally requires the use of a high 
performance spectrum analyzer with a tracking generator 
output 0] . This output provides a voltage signal of con- 
stant peak to peak amplitude over the frequency range 
of interest. Feeding this through a scries high impedance 
generates an rf current in the tank circuit which, in turn, 
creates an rf flux ( peak amplitude flip) at the SQUID 
ring. As with our previous work on SQUID ring-tank 
circuit dynamics in the inductive regime 5, 7], the de- 
tails of this behaviour depend crucially on the value of 
the static bias flux $xstat- 

In figure 3 we show four frequency responses for uni- 
directional sweeps, from low to high frequency, for the 
SQUID ring-tank circuit system of figure 2 (at T = 4.2 
Kelvin) . The data were taken at rf flux amplitudes flip 
ranging from 0.02 to 2.5$o, with <i>xstat = {n+ 1/2) $0. 
These frequency responses were recorded using a high 
performance, but conventional, commercial spectrum an- 
alyzer (a Rohde and Schwarz FSAS system) which could 
only scan in this unidirectionally manner. What is clear 
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FIG. 3: (a) to (d) T = 4.2 Kelvin, $stat = (n + 1/2) $o plots 
of the unidirectional frequency sweeps (from low to high) for 
the weakly hysteretic (/? = 2.26) ring-tank circuit system of 
figure 2 for different peak rf flux amplitudes at the SQUID 
ring (i.e. jiip — 0.02,0.5,2.3 and 2.5<&o ) showing the punch 
through at higher amplitude (e.g. jiif ^ 2<I>o). 
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from figure 3 is that as ^ip becomes a significant fraction 
of a flux quantum the lineshape of the frequency response 
curves changes dramatically. In particular, for ^ip ^ 2<i>o 
there is a punch through effect which leads to the cre- 
ation of fold bifurcation cuts. If we refer back to figure 
2, this value of peak rf flux jiip equates to a drive current 
corresponding to the feature shown bracketed in figure 2, 
but for the case of static bias flux $xstat = {n + 1/2) $0. 
The significance of this was revealed when bidirectional 
frequency sweeps (from low to high and then high to low) 
were taken using a spectrum analyzer of our own design. 
This spectrum analyzer is shown in block form in figure 4. 
As an illustration of its utility, we show in figure 5 the re- 
sult of this bidirectional sweep for the curve (d) of figure 
3 (/Lt(/3 = 2.5$o), again at $xstat — {n + 1/2) <I>o. As can 
be seen in figure 5, there are four fold bifurcation cuts, 
one each at the lowest and highest frequency excursions 
in the response, and two inner cuts, one on the sweep up 
and the other on the sweep down. Within the experimen- 
tal resolution, as reflected in the data points, these inner 
cuts appear to overlay exactly. To our knowledge, this 
has not been observed before. For this bidirectional re- 
sponse it seems reasonable to adopt the description pinch 
resonance. 

The frequency response curve of figure 5 and those 
plotted in figure 3 as a function of coupled rf amplitude, 
are perfectly typical of the low /3 value (« few) hysteretic 
SQUID ring-tank circuit systems that we have studied. 
To aid in interpretation we have provided in figure 6 a 
schematic of three different bidirectional frequency re- 
sponse curves: (a) where the inner pair of opposed fold 
bifurcation cuts are well separated in frequency (b) where 



FIG. 4; Block diagram of the spectrum analyzer used to 
record the bidirectional frequency response curves of flgure 
5. 
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FIG. 5: Expanded, bidirectional, frequency sweeps (low to 
high and then high to low) for the response (d) in the data 
of figure 3, again with jiip — 2.5$o and T = 4.2 Kelvin. This 
value of peak flip equates to a drive current, at static bias flux 
'I'xstat ~ (n -I- 1/2) $0, corresponding to the feature shown 
bracketed in figure 2. 



these cuts superimpose in frequency (i.e. in our termi- 
nology, a pinch resonance) and (c) where the fold bi- 
furcations have moved beyond one another. It seems 
clear that the data of figure 5 correspond to the bidi- 
rectional response of figure 6(b) since this response can 
obviously be differentiated from those of figures 6(a) and 
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FIG. 6: A schematic of three diflerent bidirectional frequency 
response curves; (a) where the inner pair of opposed fold bi- 
furcation cuts are well separated in frequency (b) where these 
cuts superimpose in frequency (a pinch resonance) and (c) 
where the fold bifurcations have moved beyond one another. 
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FIG. 7: Experimental bidirectional frequency response curves 
for the SQUID ring-tank circuit system of figures 3(d) and 5 
(T = 4.2Kelvin, flip = 2.5"I'o) at three values of static bias 
flux at the ring $xstat = n$o, (n -I- 1/4) $0 and (n + 1/2) $0. 



6(c). It is interesting to note that the existence of any 
one of the bidirectional response curves shown schemat- 
ically in figure 6 could not be revealed using standard 
commercial spectrum analyzers since these always sweep 
the frequency unidirectionally, i.e. from low to high. 
For completeness, we show in figure 7 the bidirectional 
frequency responses (with both up and down sweeps) 
for the ring-tank circuit system of figures 3(d) and 5 
(T ~ 4.2 Kelvin, flip = 2.5$o) at the static bias flux val- 
ues $xstat = n<^o, {n + 1/4) $0 and (n + 1/2) As is 
apparent, the pinch resonance can only be seen in one of 
these responses, i.e. at $xstat = {n + 1/2) 



III. THEORETICAL MODEL 



A. RSJ-I-C description 



In the standard dynamical description of a SQUID ring 
it is customary to invoke the well known (quasi-classical) 
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resistively shunted junction plus capacitance (RSJ+C) 
model 01 of the weak link in the ring. Here, a pair 
(Josephson) current channel is in parallel with an effec- 
tive weak link capacitance C and a normal current chan- 
nel of resistance R. For an applied current I > Ic, the 
normal channel opens up with a value of R which is char- 
acteristic of the particular weak link in the ring. In this 
(and other) models the SQUID ring moves in the space of 
total included magnetic flux $ = (E>x + A/g with an effec- 
tive mass given by the weak link capacitance (C) . Within 
the RSJ-I-C model it is always assumed that the SQUID 
ring can be treated as a quasi-classical object, such that 
C is relatively large, typically « 10^^^ to IO~^^F. 

In the RSJ+C description the SQUID ring and tank 
circuit constitute a system of coupled oscillators with 
equations of motion which can be written in the form 

tank circuit 



Rt 



Lt 



(1) 



where Rt is the tank circuit resistance on parallel reso- 
nance of the coupled system and is the drive current 
(figure I), which may contain both coherent and noise 
contributions. 
SQUID ring 



resonances, such as we have presented in figure 5, and 
which match with the experimentally determined circuit 
parameters, show up as solutions of these equations. 



1. Pinch resonances 

With a knowledge of the system parameters, and an 
estimate of Ic , we can solve the coupled equations of 
motion and numerically to find the value, at a 
given $xstat, of the rf tank circuit voltage as a func- 
tion of the frequency of the fixed amplitude sinusoidal 
drive current /in . This voltage response can then be 
used to compute the bidirectional frequency responses of 
the ring-tank circuit system. To do this we considered 
an interval [wmin/STT to ti;i„ax/27r] in the drive frequency 
which contained the (peak amplitude) resonant frequency 
of the system. This interval was divided into N bins 
which gave us a resolution of Suj = {oj max — ^min) /N in 
the frequency domain. We then drove the tank circuit at 

LJi-n (for /in)= Wniin + nSuj for n = 0, 1, , N increments 

in frequency (this being done without changing the ini- 
tial conditions of the system), thus obtaining the power 
in the tank circuit at each of these frequencies via Q 



-H/eSm( — 



M{l-K^) 
(2) 



where I ~ In sin 



27r$ 



is the phase dependent Josephson 



current for the weak link in the ring. 

Until recently little effort appears to have been directed 
to solving the coupled equations in their full non-linear 
form, understandably given the level of computational 
power required to deal with them. Instead it has often 
been assumed that it is sufficiently accurate to linearize 
the SQUID ring equation of motion. This is equivalent 
to invoking the adiabatic condition that the ring stays 
at, or very close to, the minimum in its potential 



2A 



27r 



as this varies with time. The SQUID equation Q then 
reduces to 



Ir sin 



A (I - i^2) M (I - K'^) 



(4) 



However, as we have shown in recent work |l2l |. this 
approximation should be treated with great care, par- 
ticularly in the hysteretic regime of behaviour. In order 
to avoid this problem we have chosen to solve the full 
coupled equations of motion [(Q and for the system 
using values of C and R considered typical of a weakly 
hysteretic SQUID ring. We shall demonstrate that pinch 



P (Win) 



(27rTO/aJi, 



2Trra / 



If exp (iuJi^^t) dt 



(5) 



for a given integer m. We note that m is taken so that 
the power is integrated over a few Q cycles of the tank 
circuit oscillation to reduce the effects of transients. In 
the situation where the non-linearities in the coupled 
ring-tank circuit system are strong, and where, typi- 
cally, the rf drive amplitudes are large, we have already 
shown H Qthat fold bifurcation regions (i.e. regions 
where more than one stable solution exists) can develop 
in the resonance curves of the system. 

To make a comparison with the experimental data of 
figure 5 we estimated the value of the /3-parameter from 
the 0° SQUID magnetometer characteristic of figure 2 us- 
ing the ratio (on the current axis /in) of the initial riser 
to the interval between the periodic features in this char- 
acteristic. This yielded a (3 value of 2.26, which, for a 
ring inductance of 6 x 10~^°//, set /c — I.23/i^. In a re- 
cent publication we described a quantum electrodynamic 
model |31] of a SQUID ring which has allowed us to de- 
rive an expression for the effective ring capacitance (C), 
knowing the superconducting material from which it is 
constructed and the value of the weak link critical cur- 
rent. With niobium as the material, and /c = I to a few 
PlA, this yields C ^ lO^^'^F. An ahernative model 3|, 
using a different approach, gives rise to similar size ca- 
pacitances. From the literature [13, H^, a corresponding 
effective dissipative resistance in the ring ss 10017 appears 
perfectly reasonable, setting a ring CR time constant of 
10~^2secs. 
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FIG. 8: Computed (RSJ+C model) bidirectional, frequency 
response curves (low to high and then high to low) using the 
ring and tank circuit parameters of figures 2 and 5 with Ic = 
1.23pA, C = lO'^F and R = 100^ {CR = 10"^^ sec). 

The theoretical bidirectional frequency response curves 
corresponding to the experimental data (and circuit pa- 
rameters) of figure 5 are shown in 

figure 8 , with Ic and CR taken to be 1.23/i^ and 
10~^^ sees, respectively. These were computed, sweep- 
ing forwards and backwards across the defined frequency 
window, using a 4**^ order Runge-Kutta-Merson routine 
with adaptive step size error control. The SQUID ring- 
tank circuit coupling {K^^ and Q parameters were set 
as for figures 3(d) and 5 together with jjnp ~ 2.5<f'o and 
^xstat — {n -\- 1/2) $o- From numerous computations we 
have found that the shape of the bidirectional response 
in this RSJ-I-C approach is determined principally by the 
CR product and is only weakly dependent on the particu- 
lar values of C and R chosen, at least in the ranges 10~^^ 
to 10~^^F and 1 to lOOil. With the circuit parameters 
detailed above, we found that the best fit to figure 5 (the 
pinch off resonance) is made with CR = 10^^^ sees. As 
for the experimental data displayed in figure 5, the theo- 
retical inner fold bifurcation cuts superimpose very well, 
at least to the limits of our computational accuracy. This 
is reflected in the computed data points in the region of 
these inner cuts. 

Clearly, the correspondence between the theoretically 
calculated response curves of figure 8 and the experimen- 
tal bidirectional responses of figure 5 is good, both in 
frequency range and lineshape. For example, at <&xstat = 
(n -|- 1/2) $0 the experimental and theoretical responses 
appear to meet exactly at one point in the frequency 
domain. However, we must caution that due to the nu- 
merical nature of the calculation it is not possible theo- 
retically to prove that this single frequency meeting point 
is exact. Nevertheless, by trying different amplitudes for 
/in we have found that we can narrow the separation 
between the solutions down to an arbitrarily small level. 
This reflects the situation found experimentally. We note 
that we have been unable to generate this theoretical fit 
to our experimental bidirectional reponse curves using 
the linearized SQUID equation Q). 



With regard to the alternative frequency response 
curves of figure 6, and the experimental data of figures 
5 and 7, we note that it may be appropriate to treat 
this weakly hysteretic SQUID ring-tank circuit system 
as an extremely sensitive finite state machine If so, 
we may be able to compare the transition state diagram 
for the different cycles shown in figure 6 [(a), (b) and 
(c)]. As an example, we may ask whether it is possi- 
ble in the pinch resonance (b) to go between the lower 
left and lower right branches without jumping up to the 
loop. This could happen occasionally but still be beyond 
the sensitivity limits of our present experimental appa- 
ratus. If it does occur there will be a sharp transition 
between the full pinch loop resonance [fine black line in 
figure 6(b)] and the resonance displaying sets of bifur- 
cation cuts., i.e. in the one case the resonance loop is 
not accessed, in the other (the latter) it is. In principle, 
the two processes could be distinguished since the path- 
ways around the loop involve very much larger voltages 
than the crossing between the lower branches of the reso- 
nance. Such contrasting processes would, of course, be of 
scientific interest but could also have potential in device 
applications. 



IV. CONCLUSIONS 

In this paper we have shown that a single weak link 
SQUID ring, weakly hysteretic and coupled to a simple 
parallel resonance tank circuit, can display quite remark- 
able, and previously unexpected, non-linear behaviour 
in the frequency domain. Thus, we have demonstrated 
that this system exhibits opposed fold bifurcation cuts 
in its bidirectional frequency response curves which, at 
specific values of rf and static bias flux, and within our 
experimental resolution, superimpose exactly. Since sin- 
gle fold bifurcations in the upper or lower branches of 
a strongly driven resonant system are usually taken as 
one of the standard examples of the onset of non-linear 
dynamics |^, the onset of opposed fold bifurcations is 
of special note. That this can result, as a limiting case, 
in the creation of a pinch resonance would appear to be 
exceptional and, in our opininion, of very considerable 
interest to the non-linear community. Nevertheless, as 
we have also shown, such pinch off resonances can be 
modelled well within the full, non-linear RSJ-f C descrip- 
tion of the ring-tank circuit system. Together with other 
recent published work these results would appear 
to suggest that there is still a great deal to be explored 
in the non- linear dynamics of SQUID ring-resonator sys- 
tems in the quasi-classical regime. It is also apparent 
that the new non-linear dynamical results obtained in 
this regime point to future research needed to implement 
quantum technologies based on SQUID rings (and re- 
lated weak link devices) interacting with classical probe 
circuits 0, 0, 113] . From both perspectives it is clear that 
the SQUID ring which, through the cosine term in its 
Hamiltonian, is non-linear to all orders, is set to continue 
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to reveal new and interesting non-linear phenomena of 
importance in the field of non-linear dynamics and in the 
general area of the quantum-classical interface as related 
to the measurement problem in quantum mechanics. It 
is also apparent that in both the classical and quantum 
regimes the non-linear dynamics generated in coupled cir- 
cuit systems involving SQUID rings could form the basis 
for new electronic technologies. 
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